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O I Abstract 

o 

We calculate the set of 0{a'^) corrections to the branching ratio and to the photon 
energy spectrum of the decay process B -^ Xgj originating from the interference of 
/\ • diagrams involving the electromagnetic dipole operator Of with diagrams involving 

^ . the chromomagnetic dipole operator Og. The corrections evaluated here are one of 

the elements needed to complete the calculations of the B — )• Xg^ branching ratio at 
next-to-next-to- leading order in QCD. We conclude that this set of corrections does 
not change the central value of the Standard Model prediction for Br(i? — >■ Xg-y) by 
more than 1%. 



1 Introduction 

The first estimate of tfie B — )■ Xg'-f brancliing ratio witliin tlie Standard Model at tlie next- 
to-next-to- leading order (NNLO) level was published some years ago [1]: 

Bt{B ^ X,7)sM,i<;,>i.6GeV = (3.15 ± 0.23) x 10"^ . (1.1) 

This estimate combines a number of different corrections which were calculated by several 
groups [214T2]. The prediction given in Eq. (11. ip must be compared with the current world 
averages, 



(3.55 ± 0.24 ± 0.09) x 10"^ , (HFAG) \m 
(3.50 ± 0.14 ± 0.10) X 10-^ , Pij 



Br(5^X,7)exp,£;^>1.6GeV = <! /o ^n , ni . , n . nN , , 1 n-4 n.i (1-2) 



which include measurements from CLEO, BaBar and Belle [T5HI7]. The central values of the 
theoretical prediction and of the HFAG average are compatible at the 1.2cr level, while both 
the theoretical and experimental uncertainties are very similar in size (about 7%). Since 
the experimental uncertainty is expected to decrease to 5% by the end of the B-factory era 
(which is already indicated by the average given in the second line of Eq. fll.2p ). it is also 
desirable to reduce the theoretical uncertainty accordingly. 

Unfortunately, at this level of accuracy, the theoretical uncertainty is dominated by non- 
perturbative contributions. As long as one restricts the analysis to processes mediated by 
the electromagnetic dipole operator O7 = aem/(4vr)mft {sa^^ Pjih) F^^ alone, non-perturbative 
effects are well under control [T81422] . However, as soon as operators other than O7 (such as 
the chromomagnetic dipole operator Og = 5's/(167r^)mb {sa^^'^ Pf{T°-b) G^j,) are involved, one 
encounters non-perturbative effects of O^asAQcn/'nT'b). At present, the latter can only be 
estimated [23] . Hence a 5% uncertainty related to all of the unknown non-perturbative effects 
has been included in Eq. (11. ip . A further reduction of the theoretical uncertainty below the 
5% level seems to be rather difficult |23]. Still, given the importance of Bt{B — )■ Xs'j) in 
constraining physics scenarios beyond the Standard Model [25], it is worth to reduce the 
perturbative uncertainties as much as possible. 

In particular, it would be desirable to reduce the uncertainty associated to the interpola- 
tion in rric which was employed to obtain Eq. fll.lj) [TT] . To get rid of the interpolation in rric 
in the calculation of the branching ratio is a highly challenging task and it would represent 
a clear improvement of the theoretical prediction. Indeed, considering the work that has 
been done since the publication of [1], and the work that is still in progress, an update of the 
estimate given in Eq. (11.11) will soon be warranted. Here we would like to mention that the ef- 
fects of charm and bottom quark masses on gluon lines are now completely known (provided 
that one neglects on-shell amplitudes that are proportional to the small Wilson coefficients 
of the four-quark operators O^-Oq) [261429] . Therefore this part could be removed from the 
interpolation. Also the 0(a^/3o)-effects in the (02,02), (02,07) and (O7, 08)-interference, 
which are known |30j, were not considered in [HITT]. Finally, the complete calculation of the 
(O2, 07)-interference for rric = is well underway [21]. The latter calculation in particular 
will help to fix the boundary for the rric interpolation for vanishing rric; this in turn would 
allow one to reduce the 3% uncertainty in Eq. (11.11) due to the interpolation. For complete 
up-to-date lists of needed perturbative and non-perturbative corrections to the branching 
ratio we refer the reader to the reviews 



In this paper we calculate the complete (O7, 08)-interference corrections at O(a^) to 
the photon energy spectrum dr{b — )■ X^^'^^°'^^'^'~f) / dE^ and to the total decay width r(6 — )■ 
j^partomc^^l^ ^^^^ where Eq denotes the lower cut in the photon energy. The contributions 
containing massless and massive quark loops were already presented in [T0l[30] and [28], 
respectively; the contributions which are not yet available in the literature are the ones 
proportional to the color factors Cp and CfCa- From the technical point of view, the latter 
are the most complicated to evaluate and are the main subject of the present work. 

The paper is organized as follows: In Sec. 2 we present our results for the (integrated) 
photon energy spectrum. In Sec. 3 we provide some details about the calculation of the 
corrections proportional to a'^Cp and alCpCA by analyzing the contribution of a particular 
Feynman diagram. The numerical impact of the (O7, Og) interference on the theoretical 
prediction for Bt{B — )■ ^^7) at NNLO is estimated in Sec. 4. Finally, we present our 
conclusions in Sec. 5. 



2 Results for the (integrated) photon energy spectrum 

Within the low-energy effective theory, the partonic b — > ^^7 decay rate can be written as 

r(6 ^ xr-7)i..>i.o = ^^-^[^'^^ \v,,v:f j^cf (/.)cf (/.) f\z^-^^j^, 

(2.1) 

where m^ and rnb{fi) denote the pole and the running MS mass of the b quark, respectively, 
C^^i^fi) indicates the effective Wilson coefficients at the low-energy scale, z = 2E^/mb is the 
rescaled photon energy, and Zq = 2EQ/mb is the rescaled energy cut in the photon energy 
spectrumjj 

As already anticipated in the introduction, we will focus on the function dG-/s{z, fi)/dz 
corresponding to the interference of the electro- and the chromomagnetic dipole operators 

<^7 = TTT^ ^bif^) {sa^'PRb) F,, , (2.2) 

Os = yI^ rn,{^^) {sa'^'' PrT%) G% . (2.3) 



dGjs 



In NNLO approximation G^s can be rewritten as follows, 

(2.4) 
where ^^(/i) indicates the running coupling constant in the MS scheme and 



}^ = '^C,Yi^\z,,)+('^yC,Yi^\z,,) + Oia^ 



Y^'\z,^i) 



■^ (33 - 2vr^) + 1^ In ( J^) 
9 ^ ^3 \mbj 



6{l-z) 



^In this paper we assume that the products Cj^^ {^) Cj^ (n) are real quantities. Therefore our formulas 
are not applicable to physics scenarios beyond the Standard Model which produce complex short distance 
couplings. 



,I(.^.^]-1{^-^)hi- 



z . 



(2.5) 



The function Y^'^\z,fi) can be split further into a sum of contributions proportional to 
different color factors: 

F(2) (z, /i) = C^y (2'^^) {z, /x) + C^F(2,cA) ^^^ ^) 

+ TniViF^^'"^) (;2, /x) + TniVHF^^'^^) {z,fi)+ TrNvY^'''''^ {z,fi). (2.6) 

Here, Nl, Nh and A^y denote the number of light {mq = 0), heavy (m^ = nib), and purely 
virtual (m^ = rric) quark flavors, respectively; Cp, Ca and Tr are the SU(3) color factors 
with numerical values given by 4/3, 3 and 1/2, respectively. The expressions for the functions 
y^^''-'(z,/i) with i = NL, NH, NV can be found in [28]. The main result of the present work 
are the so far unknown functions Y^'^''^\z, fj.) with i = CF, CA, which are given by 



F(2,cf)^^^^) 



64 

-37.1831 L„ 

3 ^ 



128 



Ll (5(1 -z)- 11.7874 



ln(l - z) 

l-z 



- 20.6279 



l-z 



41.7874 ln(l - z) - 6.6667 ln2(l - ;z) 



~ 64 

+ /i(2)-12F«(^,m,)L^ + yif«(^,m,)L^, 



(2.7) 



F(2'-)(^,/i) 



4.7666 + ^L, + ^Lj) 5(1-.) 

34 ~ 

- 6.5024 ln(l - z) + f2{z) + — F^^^ (^, mb)L 

o 



M ' 



where 



5 TT^ . ^,^ , 



ln(l - z) 

l-z 



l-z 



z+l ^ , 7 + z-2z^ 
ln(l — -2) H 



(2.^ 



(2.9) 



/i(^) 



20.6279 - 108.484 z + 13.264 z^ + 16.1268 z^ - 33.2188 / 
+ 69.8819 / - 111.088 / + 118.405 / - 79.6963 / + 29.929 / 
- 4.76579 z^° - 56.8265 (1 - z) ln(l - z) 



8.11265 (l-z)ln2(l 



5.77146 (l-zlln^d 



(2.10) 



f2{z) = 17.0559 z + 20.9072 z"^ - 0.471626 z^ + 10.1494 z^ 

- 17.4241 z^ + 24.7733 z^ - 20.4582 / + 8.47394 z^ - 0.173599 z^ 



0.657813 z^" + 5.66536 (1 - z) ln(l - z) 



11.1319 (1 - z) ln^(l -z) + 1.3999 (1 - z) ln^(l 



(2.11) 



Note that function H^^\z, fi) also appeared in Eq. (2.11) of Ref. |9j, and that we introduced 
the short-hand notation L^ = \n{fi/mh). 

In Eqs. (12. 7p and (12. 8 p the ^-dependence of the /^-dependent terms and of those terms 
which become singular in the limit z — t- 1 is exact. The functions fi{z) and f2{z) were 
instead obtained by making an ansatz for our numerical results of the non-singular parts, 
using the functional form 
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fi{z) = J2ci,jZ^ + c^M{l-z)Hl-z)+CiM'^-z)ln\l-z) + Q,u{l-z)\n\l-z). (2.12) 

j=0 

The coefficients cq, . . . , C13 were then determined by performing a least-square fit, using 100 
specific 'data'-points. We checked that the fit-functions remain essentially the same when 
changing the set of data-points. In particular, the integrals of the fit-functions, taken in an 
interval [^o, 1] (0 < 2:0 < 1); remain basically unchanged. The same holds true when changing 
the functional ansatz given in Eq. (I2.12p . e.g., to contain additional terms proportional to 
(1 - zf ln"(l - z), with n = 1, 2, 3. 

The plus distributions appearing in Eq. (12. 7p are defined as 



1 , rin"(l - z) 
dz 







z 



9{^) 



1-z 



[^w-^(l)] 



(2.13) 



where g{z) is an arbitrary test function which is regular at z = 1, and n = 0,1. In case the 
integration does not include the endpoint 2; = 1, we have (c < 1) 



dz 



ln"(l - z) 
1-z 



9{z) 



^ ln-(l - z) 
dz^— g[z) 



(2.14) 



We observe that the plus distributions are present only in the part of the spectrum 
proportional to C'jp, see Eq. (12. 7p . This is in agreement with the results reported in [36] : 
following the procedure presented in that work, it is possible to conclude that the plus 
distributions appearing in the (O7, Og) component of the photon energy spectrum at C(a^) 
must be the same ones as in the (O7, O7) component of the spectrum at 0{aa) (up to an 
overall factor). In particular 



y(2,cF)| 



z,mb] 



plus distrib. 



(33 - 27r2) 



ln(l - z) 



1-z 



+ 



7 



1 



-11.7874 



ln(l - z) 
1-z 



20.6279 



1-z 

1 

1-z 



(2.15) 



The structure in Eq. (I2.15P emerges in our diagrammatic calculation from delicate cancella- 
tions among several contributions, and it provides a valuable test for our result. 



3 Technical details about the calculation 

In order to obtain the results in Eqs. (12.71) and (12.81) one needs to evaluate Feynman diagrams 
contributing to the process 6 — )■ 37 up to two loops, one-loop Feynman diagrams contributing 
to the process b — )■ s'jg, and tree-level Feynman diagrams contributing to the processes 
b — >■ s'ygg (plus corresponding diagrams involving unphysical ghosts in the final state) and 
b — )■ s'yss. As discussed in detail in [9] , the interferences among the various partonic diagrams 
with 2, 3 and 4 particles in the final state are in one-to-one correspondence with the 2-, 3- 
and 4-particle cuts of the three-loop 6-quarks self-energy diagrams shown in Figs. |3l HJ [5l [6l 
and [71 provided that the cut goes through the photon propagatocl. The 2-particle cuts of the 
self-energy graphs correspond to the interference of a tree-level and a two-loop diagrams for 
the process b — )■ 57, or to the interference of two one-loop diagrams for the same process. The 
3-particle cuts split a self-energy graph into a one-loop and a tree-level diagram contributing 
to the process b — t- s'-fg. Finally, the 4-particle cuts correspond to the interference of two 
tree-level diagrams for the process b — > sgg or for the process b — )■ s'-fss. 

The contribution of each single cut to the photon energy spectrum can be obtained by 
employing the Cutkosky rules [37H39] . In order to keep the energy of the photon fixed, it is 
necessary to insert in the integrands a factor [7] 

_ Pb:pA ^ 2m,6{{p,-p,f - (1 - z)ml) , (3.1) 

nib J ^ 

where pf, and p^ denote the four momenta of the b quark and the photon, respectively. Finally, 
the delta function in Eq. (13. ip and all of the delta functions originating from the Cutkosky 
rules for the propagators crossed by a cut can be rewritten as differences of propagators as 
follows iOllill, 

1/1 1 \ 

(3.2) 



27ii \g2 _ fYi2 _ ^Q g2 _ ^2 _|_ j^q J 



It is then possible to evaluate all of the relevant cuts of the three-loop self-energy diagrams 
by first identifying a set of Master Integrals for each cut, and then by evaluating those 
Master Integrals by means of the tools and techniques usually employed in the calculation 
of multi-loop Feynman integrals. As usual, we work in c? = 4 — 2e space-time dimensions 
to regularize ultraviolet, infrared and coUinear singularities. In the rest of this section we 
discuss some features of this procedure by considering one of the simplest Feynman diagrams 
we encountered in the course of the calculation. More technical details concerning the 
parameterization of the phase-space integrals can be found in [Sl[nill2] . 

Let us consider the topology displayed in Fig. [H which corresponds to the last diagram 
shown in Fig O We carried out the reduction to Master Integrals by means of the package 
AIR ^3] and by means of a private code written by one of us (in order to have a cross check). 
In the diagram corresponding to the topology in Fig. [H only a 2- and a 3-particle cut are 
present (they are indicated by the dashed red and blue lines in Fig. [1], respectively). We will 
first concentrate on the 3-particle cut; the 2-particle cut will be discussed at the end of this 
section. Because of the delta function in Eq. (13. ip . one finds (fci = p^) 

{Pb ~ hf = ml{l - z) ■ (3.3) 



^We only display the diagrams that contribute to the functions Y'^'^'^^{z.ji) with i — CF, CA. 




Figure 1: A three-loop self-energy topology and the corresponding 2- and 3-particle cuts. 
Thin lines represent niassless propagators; the dashed line corresponds to the photon prop- 
agator; the 2- and the 3-particle cut are indicated by the dashed red and blue lines, respec- 
tively; ki, ^2 and k^ are the loop momenta. 




Figure 2: Master Integrals for 3-particle cut of the topology displayed in Fig. [II 



therefore the corresponding propagator can be immediately factored out of the integrals. 
The reduction indicates that the 3-particle cut of the topology in Fig. [T] has three Master 
Integrals, which are shown in Fig. |2l 

Those Master Integrals have only single poles in e (as long as z ^ 1), and are sufficiently 
simple to be calculated analytically by direct integration of a suitable Feynman parame- 
terization of the virtual loop and of the 3-particle phase-space. By combining the output 
of the reduction to Master Integrals with the analytic expressions of the latter, it is pos- 
sible to find an expression for the contribution of this particular cut to the photon energy 
spectrum. It is then straightforward to check that the coefficient of the single pole in e, 
the finite part, and the 0(e) term vanish in the z — )■ limit, as expected from dimensional 
reasoning [32]. On the contrary, in the 2; — )■ 1 limit this particular cut shows divergences of 
the form ln"'(l — z)/{l — z) [n = 0, 1). Such divergences give origin to the plus distribution 
functions which survive in the part of the photon energy spectrum which is proportional to 
Cp. To see this one needs to extract a factor (1 — 2;)^^ (where b is an arbitrary integer) out 
of each Master Integral, to combine it with possible factors 1/(1 — 2;) which emerge from the 
reduction procedure, and then to replace the resulting expression by using 



(l-z) 



-l+be 



be 



6{1 



+ E 



n=0 



n! 



ln"fl 



(l-z) 



(3.4) 



The relation above explains why one needs to calculate the cofactor of (1 — z)~ '^ including 
terms of 0(e): the latter terms contribute to the 0(e°) part of the coefficient of 6{1 — z). 
For the definition of the plus distributions see Eqs. (12.131) and (12.141) . 

Very often we had to deal with Master Integrals which we were not able to evaluate 







Figure 3: Diagrams with a closed gluon loop. Similar diagrams with a closed ghost loop 
are also present. The gray square indicates the operator Og, while the gray circle indicates 
the operator Of. Thick lines represent the massive b quark, thin lines the massless s quark, 
wavy lines photons and curly lines gluons. The color factor of these diagrams is proportional 
to CfCa- 




Figure 4: Same as in Fig. [3] for diagrams containing a triple gluon vertex. The color factor 
of these diagrams is proportional to CfCa- 




Figure 5: Same as in Fig. [3| for diagrams involving double gluon emission from the Og 
operator. The color factor of these diagrams is proportional to CfCa- 
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Figure 6: Same as in Fig. \^for diagrams with color factor proportional to CpiCp — Ca/'^)- 
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Figure 7: Same as in Fig. \^for diagrams proportional to the color factor Cp. 
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by direct integration of their integral representation in terms of Feynman parameters. A 
powerful tool to be used in these cases is the differential equation method [miHTlH^ HS]. 
The goal of the method is to employ the output of the reduction procedure for a given 
topology to build differential equations which are satisfied by the Master Integrals of that 
topology. In our case, we consider differential equations with respect to z, the only variable 
which appears in the Master Integrals. In order to generate these differential equations, it 
is necessary to take the derivative of the integrand of a given Master Integral with respect 
to z. It is interesting to observe that the only factor in the integrand which depends on z is 
the propagator associated to the delta function given in Eq. (13.10 . 

A weakness of the differential equation method is the fact that there is no general strategy 
which allows to fix the integration constants which are not determined by the solution of 
the differential equations. In a number of cases it was possible to fix the missing constants 
by directly evaluating the Master Integral after setting z = from the start. However, we 
also found Master Integrals which are singular in z = 0; they appear in the reduction of 
the 3- and 4-particle cuts of the 3rd and the last diagram given in 4th line in Fig. |6l To 
fix the integration constants in those cases we exploited the fact that the coefficient of each 
term in the e expansion of each single cut contributing to the photon energy spectrum must 
vanish in the z — )■ limit. However, we had to deal with cuts for which this procedure did 
not provide enough conditions to fix all of the unknown integration constants. Therefore 
we calculated some of the Master Integrals for z = 1 to reduce the number of unknown 
integration constants; subsequently we fixed the remaining ones by considering the z — )■ 
limit of the cuts involving those integrals. Another procedure which we employed to reduce 
the number of unknown integration constants consists in integrating some Master Integral 
over z from to 1. With these methods we were able to obtain analytic expressions for 
all of the poles in e appearing in the calculation of the various cuts. For a few cuts we 
calculated the finite parts only numerically. Some of the diagrams with non-trivial endpoint 
behavior were checked by an independent calculation of the integrated spectrum. We would 
like to mention that for the (O7, 08)-interference the endpoint singularities are only present 
in 3-particle cuts of the 6-quark self energies; 4-particle cuts are free of endpoint singularities 
since the 4-particle phase space is proportional to {1 — z). 

Now we will turn briefiy to the evaluation of the 2-particle cut in Fig. [T] The reduction to 
Master Integrals is carried out along the lines of the reduction of the 3-particle cut discussed 
above. However, the insertion of the delta function given in Eq. (13. ip in the integrand is 
not necessary, since the 2-particle process takes place at fixed photon energy, z = 1. We 
calculated the four Master Integrals which appear in this case by using a numerical method 
based on sector decomposition |1S]; this technique allows to disentangle overlapping infrared, 
coUinear and ultraviolet divergences. We applied this numerical method in order to evaluate 
all of the Master Integrals arising from 2-particle cuts. However, due to the presence of 
internal thresholds, the integration over some of the Feynman parameters can only be done 
numerically after a suitable contour deformation [17H15] . 

The infrared and coUinear divergences appearing in the 2-particle cuts will cancel after 
adding the 3- and 4-particle cuts, which also suffer from infrared and coUinear divergences. 
The remaining divergences are ultraviolet. They are removed by adding counterterm dia- 
grams with appropriate Z- factor insertions (see App. [B|) . 
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4 Estimating the numerical impact on Br(S — )- Xs^) 

In this section we investigate the numerical size of the (O7, 08)-interference at O(a^) at the 
level of the branching ratio of the decay process B — ?■ Xg'-f. In order to do so, we adopt the 
notation and conventions introduced in [11]. The 0{a'j.) correction to the function Kfg{Eo, fi) 
in Eqs. (2.6) and (3.1) of [TT] is given by 



Kif{Eo,fi) 



2 



dz\Y'^^\z,fi) 



C 



p—ard{l-z) 



-C, 



41 



-27r^ 



12 L, 



r«(z,/i) 



(4.1) 



with a-j = 0.22. Note that Eq. (14. ip refers to the IS'-scheme for the 6-quark mass. The value 
used for this parameter in the numerics is ml^ = 4.68 GeV. 

Combining the results of our present paper with those of p8J we are now able to write 
down the complete expression for K-^g containing all abelian and non-abelian contributions 
as well as the effects of the masses of the u, d, s, c and b quarks running in the bubbles 
inserted in gluon lines. This complete term affects the branching ratio by an amount 



ABr(5 ^ Xs-f)E-,>Eo = MB -^ X,eu] 



cxp 



VtbV* 



Vr 



cb 



Qae 



n 



c 



AP(Eo) 



where 



AP(Eo) = 2Cj°)''^(/.)Cr^^(^) 



Qs(/i) 
47r 



.(2) 



Krs'iEo,fi) 



(4.2) 



(4.3) 



and C is the so-called semileptonic phase-space factor. In order to compare with Ref. |11] . 
we employ the numerical value for C which was obtained from a fit of the measured spectrum 
of the B -)■ XcW decay in the IS schemcl [50tl5T]. 

One might think that ABr(_B — )■ Xs'y)E^>Eo in Eq. (14. 2 p simply represents the shift due 
to K^g of the theoretical prediction given in Eq. (II. ip . This is, however, not the case because 
an approximated version of Kj^ was already included in [1]: While the /3o-part of K^s (i-e. 
-^78 when following the notation of [11]) was fully taken into account, the remaining piece, 
Kjg '^, was calculated for large values of p = ml /ml and then interpolated (combined with 
contributions not related to the (O7, 08)-interference) to the physical value of p. To remove 
from the interpolation procedure and to replace it by the exact result obtained by us. 



K. 



(2)rem 
78 



is beyond the scope of the present paper; this issue will be correctly treated in a systematic 
update of Eq. (II. ID in the near future. To get nevertheless an idea of the numerical size 
of the 0{al) contribution of the (O7, 08)-interference at the level of branching ratio for 
B — )■ Xg'-f, we can ignore this issue and simply discuss a few numerical aspects of the 
quantity ABr(i? — )■ Xs'y)E^>Eoy based on i^^g which we calculated in this paper. 

Fig. [8] shows ABr(i? — )■ Xs'y)E-f>Eo as a function of zq for /i = 2.5 GeV, 0^(2.5 GeV) = 
0.271, cf^'*^(2.5GeV) = -0.369, C^"^"*^ (2.5 GeV) = -0.171, A^l = 3 and Nh = Ny = 1. 



^ Using the numerical value for C as obtained from a fit in the kinetic scheme raises the central value 
given in Eq. (II. ip by approximately 3% [S5]. 
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Figure 8: ABr(i? — )• Xs^^e^>Eq as a function of Zq. See text for more details. 



(Note that the scale fi = 2.5 GeV defines the central value of the branching ratio given in 
Eq. Ol.ip .) The remaining numerical input parameters are taken from [H]. The thick red 
solid line shows the complete result, the thin red solid line corresponds to the contribution 
proportional to Cp (including the numerical value of Cp), the green dashed line indicates 
the contribution proportional to CpCa (again including CpCa), and the blue dotted line is 
the contribution stemming from massless and massive quark bubbles. The dash-dotted blue 
line indicates the corrections obtained by applying the large- /3o approximation |53l|5l]. The 
vertical line indicates the value of zq corresponding to the choice Eq = 1.6 GeV. It can be 
seen from the figure that the contributions proportional to Cp and CpCa cancel each other 
over almost the whole range of 2:05 resulting in a contribution of size similar to the one due 
to the fermionic corrections. Only in the region very close to the endpoint the contribution 
proportional to Cp dominates, due to its singular behavior for zq ^ 1. We stress, however, 
that the cancellations mentioned above refer to a value of /x = 2.5 GeV; they do not occur 
anymore when going to smaller values of fi. 

For a photon energy cut-off Eq = 1.6 GeV, as the one employed in Eq. 01. ip . we find the 
following numerical value for the quantity in Eq. 04. 2 p (using /i = 2.5 GeV): 



ABr(5 ^ X,7)s^> 



C, 



1.6 Gcv = — (52.21 Cp - 23.57 Ca - 2.66 C| - 2.28 Nl 

+ 8.73 Nh - 1.42 Nv) x 10"^ = -3.57 • 10 



(4.4) 



where in the last step we inserted numerical values for the color factors and we set N^ = 3 
and Nh = Ny = 1. By comparing this with the central value of the estimate given in 

Eq. (II. ip . one sees that the 0{a'^) corrections i^fg have an impact of —0.11% at the level of 

(2) 
the branching ratio. For the analogous effect generated by the large-/3o approximation of K^g 

only, Kjg , we find +0.56%. These results are, however, very strongly dependent on the 

scale fi. For fi = 1.25, 2.34, 5 GeV the corresponding numbers read —5.15, —0.23, —0.07% 

(full) and +0.10, +0.49, +1.15% (large-/3o approximation), respectively. From these results 

we conclude that the large-/3o approximation does not provide a good estimate of the full 

0(af) correction of the (O7, 08)-contribution for fi G [1.25, 5] GeV. As already mentioned, a 
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(2) 

more detailed analysis of the effect of the complete calculation of K)^ on the central value 
of Eq. (11. ip would require to repeat the interpolation procedure of [11] . While this is beyond 
the scope of the present work, we can conclude that the correction originating from the 
(O7, 08)-interference at 0{a1) will not alter the central value of Eq. (II. ip by more than 1%. 



5 Summary and conclusions 

In the present work we calculated the set of the 0{a^^ corrections to the partonic decay 
process h — )■ Xs7 which originates from the interference of diagrams involving the electro- 
magnetic dipole operator O-j with diagrams involving the chromomagnetic dipole operator 
Os. These corrections are one of the elements needed in order to complete the calculation 
of the branching ratio for the radiative decay B — )■ X57 up to NNLO in QCD. 

To carry out the calculation, we mapped the interference of diagrams contributing to 
the processes h — )■ 57, h — )■ 575', h — )■ s-^gg and h — >■ 5755 onto 2-, 3- , and 4-particle cuts 
for the three-loop 6-quark self-energy diagrams which include insertions of the operators O-j 
and 0%. Subsequently, we evaluated each single cut by employing the Cutkosky rules. From 
the technical point of view, the calculation was made possible by the use of the Laporta 
Algorithm [55] to identify the needed Master Integrals, and of the differential equation 
method and sector decomposition method to solve the Master Integrals. 

From the phenomenological point of view, it is interesting to estimate the effect of these 
corrections on the theoretical prediction for the B — )■ Xg'-y branching ratio. Our conclusion 
is that they will not change its central value given in Eq. (11. ip by more than 1%. 

At present, the largest theoretical uncertainty affecting the prediction in Eq. (II. ip is 
of non-perturbative origin. It is expected to set a lower limit of about 5% on the total 
theoretical uncertainty for the prediction of the B — )■ Xs'^ branching ratio in the near future. 
The perturbative 0{a1) corrections of the (O7, 08)-interference presented in this paper are 
a further contribution to make the perturbative uncertainty negligible with respect to the 
non-perturbative one. 
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A Two-particle cuts 

To obtain the contributions to the functions Y^^', Y^^''-^'^' and y(2,CA) (^^^^ -g^g^ (\2A\j and 
(12.61) ) originating from the 2-particle cuts of the 6-quark self energies discussed in the paper, 
we first calculated the ultraviolet renormalized on-shell matrix elements 



(O,) ^ (37IO.I6) 



7, 



For the (O7, 08)-interference, (O7) and (Og) are needed to one-loop and two-loop accuracy, 
respectively. By considering only the terms proportional to Cp and CfCa to O(a^), which 
are the ones of interest in order to obtain Eq. (12.71) and (12.81) . we write 



m = (o,), 



S. + I^C, D« + (^) C, (C. A? + CaDS + ...)+ Oial) 



(A.l) 



Note that the operator O7 in (07)tree contains the 6-quark running mass rribdj). For O7 the 
relevant results are 



D- 



(1) 



1 



(2 L^ + 2.5) -211-7 Lf,- 6.8225 



e (1.3333 Ll + 7Ll + 13.6449 L^ + 13.4779) 

e^ (0.6667 LJ + 4.6667 LJ + 13.6449 LJ + 26.9559 L^ + 26.1412) 



(A.2) 



and for Os one finds 



D 



(1) 



2.6667 L^ + 1.4734 + 2.0944 i 

+ e [2.6667Lj + 2.9468L^- 1.1947+ i (4.1888 L^ + 4.1888 
+ e^ [1.7778 Ll + 2.9468 L^ - 2.3894 L^ - 5.5373 
+ i (4.1888L2 + 8.3776 L^ + 2.1627)], 



D 



(2) 
8F 



D 



(1) 



11. 

----(2L^ + 2.5) 



5.3333 L-'^ - 32.2802 L; - 50.9612 L^ 



D 



(2) 
8A 



-1.8875 - i (4.1888 LJ + 31.4159 L^ + 29.8299) , 
15.111 Ll + 31.6617 L^ + 2.38332 + i (23.7365 L^ + 28.0745) . 



(A.3) 



Taking into account the phase-space factors in 0? = 4 — 2e dimensions, we easily obtain the 
contributions to the functions Y^^\ y(2,cF) ^^^^^ y(2,cA) -yyj^JQ}^ originate from the 2-particle 
cuts only: 



yW 



Y 



2-cutsV^' f^) 
(2,Cf), 



2 , ,, 16 ■ 

- (33 - 2n') +-L, 



6{l-z) 



2— cuts 



;2:,/^) 



-^ (-5.89368 - 10.6667 L^) + - (-15.849 - 72.6954 L^ - 53.3333 LJ) 
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+ 3.01591 - 246.244 L^ - 335.42 Lj - 135.111 Lj 6{1 - z) , 

Fa-cute (-2, Ai) = [4.76664 + 63.3235 L^ + 30.2222 LJ] 5(1 - z) . (A.4) 

The results given in this appendix were aheady used by one of us in [57] . 



B Renormalization constants 

In this appendix, we collect the explicit expressions of the renormalization constants needed 
for the ultraviolet renormalization in our calculation. The strong coupling constant is renor- 
malized in the MS scheme: 

«r^=/i'^(|^)'^r«.(/x), (B.i) 



with 



7MS_i ^ (^^ n 4^ ,^ \ a,(/i) 2N 



and Np = N^ + N^ + Ny- The 6-quark mass which appears in the operators O-jg, as well 
as the Wilson coefficients Cy^ themselves, are also renormalized in the MS scheme |58j : 



zW=l + li^Cp- Ca) ^ + 0{al) . (B.3) 

All the remaining fields and parameters are renormalized in the on-shell scheme. The on-shell 
renormalization constant for the 6-quark mass is given by 

The renormalization constants for the gluon field and the s- and 6-quark fields are 

2e 
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Z,T = l-C^r(e)e-^fii)"^ + OK^), (B.5) 

1 — ze \mbj Att 

where p = ml /ml. 
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